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Abstract 

Q _ We confirm and expand on work by Chafin and Schafer on hydrodynamic fluctuations in the 

00- 

unitary Fermi gas. Using the result for the equation of state from an recent MIT experiment, we 
derive lower bounds for r]/n and r]/ s as a function of temperature. Re-analyzing recent quantum 
' Monte Carlo data for the shear-viscosity spectral function we point out a possible resolution for the 



tension between the viscosity bound rj/n ^ 0.3 from Chafin and Schafer and the quantum Monte 
Carlo results ry/n ^ 0.2 from Wlazlowski et al. near the critical temperature. 
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I. MOTIVATION 



Hydrodynamics allows for long lived shear and sound waves. Fluctuations can populate 
these modes and interaction of these waves can dissipate momentum. This contributes to 
the shear viscosity and other transport coefficients. 

Our study builds upon results found in the context of relativistic fluids, specifically Ref. ^ 
which was recently applied to the case of the Fermi gas at unitarity by Chafin and Schafer 
in Ref. In Ref. j2| it was found that the dimensionless combination shear viscosity over 
density rj/n must be larger or equal to 0.32 for the unitary Fermi gas at the superfiuid- 
normal phase transition temperature T = T^. Thus, there is a tension between this result 
and the result for rj/n found in the quantum Monte Carlo simulations by Wlazlowski et al. 
in Ref. jsl, who find Ti/n\j,^j,^ ^ 0.2. This provides the motivation for us to re-derive and 
generalize the result by Chafin and Schafer as well as re-analyze the results by Wlazlowski 
et al. 

II. DERIVATION OF THE CORRELATOR 

We will be interested in a particular retarded correlator of the energy-stress tensor T^'^ , 
namely 

GT\x-y) = {T^y{x)T^y{y))n, 

because the Fourier transform of this correlator in classical hydrodynamics (meaning in the 
absence of fiuctuations) can be related to the value of transport coefficients. Specifically, for 
the case of a conformal system with vanishing bulk viscosity coefficient, one has ( 2I, cf. j4|) 

GT\i^. k) = P - zr/a; + t^tru' - ^k^ + 0{uj\ k^) , (1) 

with 7] the shear viscosity coefficient and r^? the relaxation time. It is possible to derive 
this result within the stochastic formulation of hydrodynamics (see e.g. Ref. sl for a recent 
review). The key idea is to supplement the energy-stress tensor with a noise term ^^'^ that 
is modeling fiuctuations, 

T^" ^ T^i + r , {r{x)C^y)) « S\x - y) , 

where the index "cl" refers to purely classical hydrodynamic expressions and () here denotes 
the noise average (see appendix |A] for an illustration). In the case of the G^^^^ correlator, it 



has been shown in Ref. [l| that to leading order in fluctuations, the stochastic formulation 
reproduces the —iurj term in Eq. ([1]), while to next-to-leading (one-loop) order, the result 
can be written as 

duj' f d^p 

' X 



(27r) 
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A-(u;', p)A^^(a; - uj\ -p) + A^^(a;', p)A^/(u; - u\ -p) + 
A-(a;', p)Af (a; - uj\ -p) + A^^(a;', p)Af (c^ - c^', -p)) , (2) 

where r symmetric and retarded version of the velocity- velocity correlator (see 

appendix In the zero momentum limit |k| — )■ 0, the integrations in the one-loop contri- 
bution to the correlator are straightforward (cf. l|) and one finds 

IGT't"" « Pn,.,. |k| = 0) = -v^'^^ + (1 + + ■ ■ ■ . (3) 

where 7 = ^ and Pmax is the cut-off for the effective theory of hydrodynamics. Note that 
terms not included in Eq. ([3]) are those that are suppressed by higher powers of Pmax7, which 
is the "small parameter" of the effective theory. Similarly, higher loop contributions to G^^^ 
are suppressed by higher powers of Pmax7- Thus, up to first order in powers of Pmax7, the 
result for the retarded correlator is given by the sum of the "tree-level" result Eq. ([1]) and 
the one- loop result Eq. (jS]). 

III. MINIMUM ON VISCOSITY 

Using the Kubo relation 

, -lmGT\^, k = 0) 
7] = lim , 

uj-i-O CO 

the result for the physical viscosity including fluctuations to first loop order is given by 

120?7ci7r^ 

where again the index "cl" refers to "classical" contributions in the sense of turning off all 
fluctuations. As was pointed out, our calculation is systematic if 

Pmax ) 
Vcl 



so we parametrize pmax by 



Pmax = — $ 
Vcl 



where $ is a dimensionless function of thermodynamic variables that should be small. For 
instance, we may take^ $ to be a function of pressure over density, $ = $(P/p). Then the 
result for the physical viscosity may be cast into a result for the dimensionless viscosity over 
density ratio, 

V\ _ fV\ ^ 17Tp2$(p/p) 



This expression can be extremized, finding that the physical viscosity over density expression 
may never decrease below the minimum value 

n) - U/min V 1607r%3 I ' ^ ' 

The minimum value for viscosity is dimensionless and presumably universal, meaning that 



for a unitary Fermi gas it should only depend on the combination T jTp with Tp 
Assuming a power-like behavior, we have 

/l53Tp2$(p/p)\ (T 



(37r2n)2/3 



2m 



V leoTT^n^ ; \Tp, 

with a a constant. Similarly, we take $ to be given by a power law behavior with 

$(P/p) oc {Pipf , 

with /3 a constant. As shown in appendix [HI the pressure is given by the relation P 
nTG{T /Tp) and ImiT^ooG = 1. With this, we find 

n \m/ \ip/ 
which leads to the requirement a = l + /3,2 — (3 = a and hence 

2' ' 2 

Thus we obtain for the final form of the lower viscosity bound 

1/3 



^ In principle, $ may also depend on the combination T/Tp with Tp the Fermi temperature (see below). 
However, we will drop this dependence in favor of a constant that may vary around unity. 
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with k an arbitrary constant of order unity. One finds 

1/3 



nJ \n 



( 459^ \ 
^ 23/2160 

V / 



which matches the result found in Ref. 2|. 



IV. WIDTH OF SPECTRAL FUNCTION PEAK 

Collecting all terms of the retarded correlator G^j^^^{uj) from Eq.([T]) and Eq. ([2]), the full 
result for the one-loop spectral function ri{ijj) to leading order in frequency is given by 



r]{uj) fri\ 51 fvV^fT\^^^ jP 



n \n/c\ l-^i^llK) \n/c\ 



23/2120 Vn/ci \TfJ V riT 



37r(28 + 3y6) /r7\-3/2 / T 



-)J;^ + 0(c.^), (7) 



23/2980 Vn/ci \TfJ\Tf 

with k an unknown constant of order unity. Let us now assume we have arrived at the 
viscosity coefficient lim(^_!.o by some other means, e.g. a first principle lattice calculation. 

The viscosity coefficient then fixes the value of rjd/n and hence Eq. ([7]) gives the low- 
frequency behavior of the spectral function. We may then compare this result to the spectral 
function obtained by the first-principles method. 



Let perform this analysis on the spectral function for t 



le shear viscosity obtained with the 



In this article, the authors ana- 



Path Integral Monte Carlo (PIMC) approach from Ref. 
lyzed data using the Maximum Entropy Method (among others) based on the 5-parameter 
model class 

C 

r]{u, m, a, c, ai, ^2) = f{uj, «!, «2)— — 7^ + (1 - /(^, "1, "2)) N{uj, m, a, c) , 

where / = e-°^°^ ^^Xi^-is) ^ ^ = -^e-'-"^^/'^/''\ to = Lo/Tp and C is Tan's contact 
(which is determined by other means and hence will be assumed to be known). This class 
of models implements the known large frequency behavior because f{u) 00) — i- 1. Inter- 
estingly, note that the model also is compatible with the low frequency behavior in Eq. 
because f{u ^ 1) = It is therefore possible that this model reproduces the cor- 

rect low-frequency behavior, even though it is not guaranteed. We re-analyze the data from 
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a 



Ref. [3| for the correlator G{t) for the = 8, A^t = 52 simulation on T/Tp = 0.355 using 
a simple minimization of 

r^^ ^ 1 .^cosh{u{TTF-^)) ^ 

G[t)-- dujr]{u,m,a,c,ai,a2)uj — . ^0, (8) 

7r JO smh[u)I f/ 2/ 1 ) 

neglecting the first three values of r for reasons explained in Ref. |3|]. We find a very small 
X^/dof ~ 4 X 10~^, but also observe that the minimum is very broad in parameter space 
{x^ does not change appreciably when changing parameters by order one). The extracted 
spectral function is shown in Fig. [T]^a). As can be seen from this figure, in our re-analysis 
of the PIMC data, the extracted shape of ri{oo) is close to the result reported in Ref. j^, 
but the tail asymptotics are only matched at much higher values of frequency. Nevertheless, 
the low frequency behavior, including the extracted viscosity coefficient r] = \im^^Qri{u) 
is rather similar. Using the viscosity coefficient to extract a value for quantity r]c\/n as 
described above (assuming ~ 1), we can then compare the PIMC results to the purely 
hydrodynamic result Eq. (jTj)^. We find that for the case of T/Tp = 0.355, the spectral 
function shape extracted from the PIMC data agrees reasonably well with the behavior 
demanded from Eq. ([7]). Because it matches the known low frequency behavior, we conclude 
that this provides additional evidence that the extracted value of the shear viscosity of a 
unitary Fermi gas at T/Tp = 0.35 is r^/n ~ 0.7. 

Let us now repeat our analysis for the low temperature case T/Tp = 0.15. Simple 
minimization of Eq. ([8]) leads to the result shown in Fig. [T](b). Comparison to the result 
from Ref. 3|] now indicates that both analysis agree with each other for large frequencies, 
while there is considerable discrepancy at low frequencies. While Ref. ^] finds a spectral 
function that is essentially featureless at low frequency, our re-analysis of the same data 
show a peak at u; = 0. Trying to convert the viscosity coefficient ri{u = 0) into a value for 
T]c\/n we find that in both cases, the constant k has to be less than unity. In particular, 
for our re- analysis we need k ~ 0.5, while to reach the viscosity value found in Ref. jsl, 
k ~ 0.05 is needed. The latter value is a much larger deviation from unity than could 
be expected. Nevertheless, with these choices for k, we can evaluate the low frequency 
behavior of the spectral function, finding the results shown in Fig. IH^b). For our re-analysis. 



^ In general, the extraction of ijci/n will not be single- valued. While negative values of rjci/n are clearly 
disallowed, we choose the smallest non- negative value of r}c\/n. 
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FIG. 1: Spectral function r]{Lo)/n extracted from PIMC simulation Ref. [3| with Nx = 8, Nr = 52 
for T iTp = 0.355 (plot (a)) and T jTp = 0.15 (plot (b)). Dots: (approximate) result published in 
Ref. [3.]. Dashed: our re-analysis of the same data using a simple fit (Eq. dE]))- Also shown is 
the known UV behavior ("high w asymptotics"). Solid lines are low frequency behavior extracted 
from the known hydrodynamic behavior Eq. ([7]). Since the T/Tp = 0.15 low frequency peak is 
extremely narrow, we have highlighted it using an ellipse in plot (b). 

the low frequency shape is roughly consistent with the hydrodynamic behavior, although 
the agreement is somewhat less than satisfactory. 

For the results for the viscosity coefficient from Ref. j^, we find that the hydrodynamic 
behavior in Eq. (I7l) results in an extremely narrow peak, in complete disagreement with the 
spectral function shape found in Ref. [s]. We suspect that if indeed such a narrow peak is 
present in the physical spectral function, it could have been missed by the method used in 
Ref. Sj to reconstruct the spectral function from the imaginary time correlation function. If 
we simply add on such a low frequency peak to the broad spectral function found in Ref. jsl, 
we find that the resulting value for the shear viscosity at T jTp = 0.15 is T]/n ~ 0.35, 
consistent with our re-analysis of the PIMC data, and consistent with the viscosity bound 
found in Ref. li. 
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FIG. 2: Summary of results for r]/n (left) and rj/s (right). The viscosity limit from Eq. ([6]) was 
used to delineate a "disallowed region" of viscosity values using the MIT equation of state from 
Ref. for P/n/T and S/N. Also shown are results from kinetic theory at high temperatures and 
the phonon scattering at low temperatures (see text for details), as well as results from Enss et al., 
Ref. Q]. Furthermore, we indicate the results from the PIMC calculation by Wlazlowski et al. Q] 
and our own re-analysis of the same data. Note that the resulting value for r]/n is consistent with 
the "upper bound" reported in the supplement to Ref. 

V. DISCUSSION AND CONCLUSION 

Employing the framework of hydrodynamics in the presence of fluctuations, the low 
frequency behaviour to the shear viscosity spectral function has been calculated. The calcu- 
lation gives a bound on the physical shear viscosity, Eq. ([6]), which matches the result found 
in Ref. ^ but disagrees with the low-temperature results from Ref. jsj . We pointed out that 
the reason for the discrepancy can be tracked to the low-frequency shape of the spectral 
function extracted in Ref. jsl. Specifically, assuming that the shape from Ref. jst is correct 
for u > O.lTp, hydrodynamic fluctuations predict a narrow but high peak at u < O.lTp, 
which effectively doubles the viscosity value. The situation is summarized in Fig. [21 where 



the viscosity bound ([6]) is shown together with the viscosity in kinetic theory (cf. Ref {g]) 



457r^/2 / 7- \ 3/2 . J. s 3/2 
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as well as the viscosity of plionon scattering in the superfluid phase, (cf. Ref 



0; 



phonons 



5.3 X 10-^/^ f^)' ^ 1.76 X 10-« (^V . (10) 



We conclude that once the low frequency peak in the spectral function is included in the 
analysis, the viscosity value from the PIMC method no longer is in conflict with the lower 
viscosity limit (Q. Also, note that in this case the viscosity value are then much closer to 
the result from Enss et al. reported in Ref. 



A 
Ref. 



so shown in Fig. [2] are results for the minimum rj/s, using the experimental data from 



for the entropy per particle as a function of temperature. The corresponding results 
for kinetic theory and phonons are 



2.768 I ^^"^ 



KT 



Tf) | + ln[3^(T/T^)3/2 



s 



phonons 



5.94xl0-«(^) . (11) 

One observes that the viscosity result from kinetic theory enters the disallowed region at 
around T/Tp ^ 0.3. We interpret this as meaning that kinetic theory looses its validity at or 
above this temperature. Similarly, the viscosity result from the superfluid analysis enters the 
disallowed region for T/Tp ^ 0.1. Again, we interpret this as meaning that the calculation 
looses validity at or below this temperature. The results for the shear viscosity by Enss et 
al. stay within the allowed region for all temperatures shown, but it is interesting to note 
that the rj/s values at low temperatures stop at the boundary to the disallowed region. 

We conclude that the theory of hydrodynamic fluctuations can provide a relevant tool to 
study the properties of unitary Fermi gases. Specifically, we expect that our findings may 
help in the analysis of future quantum Monte Carlo calculations. We plan to further develop 
this theory in the future. 
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Appendix A: Finding the Velocity- Velocity Correlator 



Let us start by writing down the Navier-Stokes and continuity equations 

^W+g^=0, -p+_(p.,)=0, (Al) 

where Iljyfc is the standard viscous fluid stress tensor. For consistency with conformal fleld 
theories the bulk viscosity, ^, has been set to zero. This is a good approximation of cold 
atom systems. Ultimately, we would like to quantify how thermal fluctuations dissipate 
momentum, acting to increase viscosity. This requires flnding the velocity-velocity correla- 
tor, so that the retarded correlator of two stress tensors can be found. In order to flnd the 
velocity- velocity correlator, a stochastic formulation of fluid dynamics is needed. To do this 
a noise term, Sjj, is added to the stress tensor, flj^ — )■ Hj^ + Sj^. The choice of noise term 
will dictate what form of fluid dynamics is found. To reproduce the Navier-Stokes equations, 
the noise must have a short distance Gaussian correlation. 

= (A2) 

(eT") = r]T{bahra + hAm)K^ " x ) (A3) 

Now with this stochastic formulation we can flnd the hydrodynamic correlation functions. 
Consider a system in equilibrium, allow for small perturbations away from equilibrium in v 
and p — >■ po + 5p. The equation of fluid mechanics to flrst order in the perturbations are 

d d Ov^ 

-^{pvk) + -^{nik + Eik)^0, dtp + po-g^ = 0. (A4) 

The velocity can be decomposed into longitudinal and transverse parts satisfying 

V-Vt{r,t)^0, V xvi{r,t) ^0. (A5) 

These restrictions give a new system of fluid dynamic equations. The linear form of the 
continuity equation yields 

^ + p,V.vi^O (A6) 
and momentum conservation yields the following system of equations 

—Vt = y^Avt , po-Q^vi = -V5P + -r)V{V ■ vi) . (A7) 
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Assuming that the velocity v{r, t) goes rapidly to zero at large distances, we can define 
the Laplace-Fourier transform 

v{k, u)^ jdr duj e-'^-'+''^%{r, t) (A8) 

In an infinite medium the Laplace-Fourier transform allows us to solve the equations of 
motion. Considering only the transverse velocity 



—vjj -\- jk'^ 

one has 



with 7 = Tj/p, the kinematic viscosity. It follows that 

<"H>-^(%-t)^. (All) 

Repeating this process for the longitudinal velocity, the entire velocity- velocity correlator 
is found to be 



2T 

{ViVj) = — 



A^- (A12) 



where 7 = §7 and we have introduced the symbol A5 for the symmetric velocity-velocity 
correlator. Now let us calculate the retarded velocity-velocity correlator. The classical 
Fluctuation-dissipation theorem reads 

X"(uj) = ^^ujC(u) (A13) 

with (3 — 1/T, being the inverse temperature and 

Qj{t) = {Ai{t)Aj{Q)), = {dAi{t)dA,{0)) . (A14) 

This is the Kubo or relaxation function with the subscript c stands for cumulant or con- 
nected. In Fourier space Cij{uj) is the symmetric correlator of the fluctuations. 

To find the retarded correlator, we must examine the generalized susceptibility. It is 
related to the time derivative of the Kubo function by 

Xij{t) = -(3e{t)Cij{t) (A15) 
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FIG. 3: Left: results for the normalized pressure P/Po from the MIT experiment Ref. [9|] and a 5 
parameter fit (see text for details). Right: The function G = as a function of T/Tp. 



where Q{t) is a step-function function. The generahzed susceptibihty x can be identified 
as the retarded correlator. Utihzing the relation x" = Ini(x) and that x is analytic in the 
upper half plane we can find the retarded velocity correlation function: 



(A16) 



The Kubo-Martin-Schwinger conditions relate the symmetric and retarded correlator in 
quantum mechanical regimes. In Fermi systems the response functions are related to bosonic 
operator (e.g. bilinears of fermi operators). Therefore, a classical version of the KMS relation 



is found to be G. 



2T 



Appendix B: Equation of State 

The equation of state of a cold unitary Fermi gas was measured experimentally by the 
MIT group [9]. The result from Ref. joj] for the pressure (versus fugacity ( = e~^^^) is 
shown in Fig. [3t^a). The result for the pressure is normalized by the ideal non- interacting 
single-component Fermi gas pressure, 

Po(/i,T) = -rA,iLi5/2(-r'), 
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where X^b = y ^ is the deBroghe wavelength. Also shown in Fig. [3]^a) is a fit of the 
experimental data using the parametrization 



MC) = P(/i,T)/Po(/i,r) 



which is a variant of the approach followed in Ref . 12| . A least square fit gives 



ci = 101.591, C2 = 158.375, 03 = 104.884, C4 = 4.30128 , 05 = 67.3894. 

Note that C4 is related to the Bertsch parameter ^ as C4 = ^"^''^ from the relation P(/i, T = 
0) = ^~^/^Po(/^, T = 0). The fitted value of C4 therefore corresponds to ^ = 0.378. Using the 
thermodynamic relation n = ^ and the identity n = A^g "^^^ '^^s 

and hence a (numerical) result for ( = ({T/Tp) by inverting the above functional relation. 
Using the above relations the pressure may be written as 

1 



P(/x,T)=nTG(C), G(C) 



-c 



For convenience, the function G is plotted in Fig. |3]^b) as a function of T/Tp- As can be 
seen, the form of G implies that the pressure is close to nT, but not monotonically larger 
or smaller than nT in the domain considered. We have checked that limy^oo G{T /Tp) — )■ 1, 
as expected. 
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